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Abstract 

Type I X-ray bursts are thermonuclear explosions that occur in the envelopes of accreting 
neutron stars. Detailed observations of these phenomena have prompted numerous studies in the- 
oretical astrophysics and experimental nuclear physics since their discovery over 35 years ago. In 
this review, we begin by discussing key observational features of these phenomena that may be 
sensitive to the particular patterns of nucleosynthesis from the associated thermonuclear burning. 
We then summarize efforts to model type I X-ray bursts, with emphasis on determining the nu- 
clear physics processes involved throughout these bursts. We discuss and evaluate limitations in 
the models, particularly with regard to key uncertainties in the nuclear physics input. Finally, 
we examine recent, relevant experimental measurements and outline future prospects to improve 
our understanding of these unique environments from observational, theoretical and experimental 
perspectives. 
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1 Introduction 

Type I X-ray bursts are the most frequent types of thermonuclear stellar explosions in the Galaxy. As 
well, they rank third after supernovae and classical novae in terms of total energy output in a single 
event, for phenomena observed in the Galaxy. Ninety-six Galactic low-mass X-ray binary systems 
exhibiting such bursting behaviour have been discovered to date^] 

X-ray bursts were discovered in 1975 [l] through observations made with the Astronomische Ned- 
erlandse Satelliet (2] of a previously known X-ray source, 4U 1820-3C(^J While a similar event had been 
observed in 1969 from the source Cen X-4 using the Vela 5B satellite [3j[4], the authors did not relate 
the observed features to a new, distinct type of emission. Consequently, that source was not associated 
with what we now call type I X-ray bursts until 1976 (5||6]. At a distance of ~1 kpc |!6], Cen X-4 is 
the nearest-known X-ray bursting source [7] and has yielded the brightest X-ray burst ever recorded 
(~60 Crab^] |6|). These initial, pioneering discoveries were soon followed by the identification of more 
bursters. While most X-ray satellites have detected X-ray bursts, great advances in the field were made 
through observations by EXOSAT (operating from 1983 to 1986 |§|), BeppoSAX (operatiing from 1996 



1 See http://www.sron.nl/~jeanz/bursterlist.html for a list of known Galactic bursting sources. 

2 Most X-ray sources are named using letters from the satellites that discovered them (e.g., 4U stands for the fourth 
catalogue of the satellite Uhuru, the first X-ray observatory) and numbers corresponding to their coordinates in right 
ascension (1820 stands for 18h20min) and declination (-30 deg) in the sky. These sources may also be named after 
the constellation where they are located and the order of discovery: 4U 1820-30 is also known as Sgr X-4, the fourth 
X-ray source discovered in the constellation Sagittarius. To complicate matters further, even within one of these naming 
conventions a single source may have two names. For example, the source 4U 1820-30 has been referred to in the literature 
as 3U 1820-30 as well. 

3 1 Crab = 2.4 xl0~ 8 erg s" 
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to 2003 (9)), and the Rossi X-ray Timing Explorer (RXTE, operating from 1995 to 2012 (To)). RXTE, 
with its large collecting area and capability for high temporal resolution X-ray analysis, has observed 



more than 1000 bursts from about 50 sources 11 



Soon after the discovery of X-ray bursting sources, it was proposed that the explosions could be 



powered by thermonuclear runaways on the surface of accreting neutron stars 12 , 13 (see Section 



3.1). The discovery of the Rapid Burster in 1976 [14], with burst recurrence times as short as 10 s, 
complicated matters. The quick succession of flashes observed from this source seemed incompatible 
with the proposed mechanism, and moreover, did not match the general pattern shown by the other 



bursting sources. A classification of type I and type II bursts was therefore established 15 , the former 
associated with thermonuclear flashes and the latter linked to accretion instabilities. The Rapid Burster 
is one of only two known sources showing type II bursts, together with the Bursting Pulsar (GRO 1744- 
28); furthermore, the Rapid Burster is the only source giving both type I and type II X-ray bursts [l~6~] . 
Type I bursts display a spectral softening during the burst tail, reflecting a decrease in the effective 
temperature and thus the cooling of the neutron star atmosphere. Typically, light curve rise times 
range between ^1 and ~10 seconds and decay times range between ^10 seconds and several minutes. 
Recurrence times may vary between roughly one and several hours. In contrast, type II bursts exhibit a 
constant temperature for the duration of the burst, which can range from ^2 s to ~700 s, with intervals 
between bursts from ^7 s to ~1 hr. 

More recently, other varieties of X-ray bursts of longer duration have been discovered. Superbursts, 
first identified from the source 4U 1735-444 by Cornelisse et al. (2000) (17, 18 , have been observed in 



about 10% of all bursters. They are characterized by light curves with exponential-like decay times of 
several hours and include an extreme case, KS 1731-260, that lasted for about 12 hours [19]. Roughly 
1000 times more energetic than a typical X-ray burst (superbursts release ~ 10 42 erg per burst), super- 
bursts have much longer recurrence times of the order of a year. The mechanism for these phenomena is 



thought to involve the unstable burning of carbon 20 22 . As well, bursts intermediate in both energy 



(releasing w 10 41 erg) and duration (^ 30 min) to typical type I X-ray bursts and superbursts have 



been observed 23 , 24 . These are thought to result from thermonuclear ignition within a layer of pure 



helium 25 , although other explanations may also be viable 26 . We will not discuss these bursts of 
longer duration or type II bursts further. 

X-ray bursts (XRBs) occur on the surface of a neutron star belonging to a low-mass X-ray binary 
(LMXB). LMXBs are binary systems where a faint, low mass (M < 1 M ) main sequence or red giant 
star transfers material onto a neutron star or black hole. The transferred matter is typically enriched in 
H, He or both H and He. Furthermore, these mass transfer episodes are driven by Roche-lobe overflows, 
leading to the buildup of an accretion disk that surrounds the compact object. (In contrast to LMXBs, 
high-mass X-ray binaries consist of a compact object capturing a fraction of the stellar wind from a 
high-mass companion star.) The maximum mass accretion rate is assumed to be set by the Eddington 
limitQ which for accretion of H-rich material onto a 1.4 M neutron star is MecLcL ~2x 10~ 8 M yr -1 . 
XRB sources often have short orbital periods, in the range of 0.2 - 15 hr; notable exceptions include 
Cir X-l, GX 13+1 and Cyg X-2, with P orb = 398.4 hr, 592.8 hr and 236.2 hr, respectively. 

As accretion proceeds onto the surface of the neutron star, the transferred material is gradually 
compressed to large densities and becomes degenerate. This compression also heats the envelope, 
creating conditions favourable to the ignition of the fuel through nuclear reactions. Since the envelope is 
initially degenerate, it cannot expand to accommodate the energy released through the nuclear processes. 
As a result, both the temperature and nuclear energy generation rate in the envelope increase. This 
drives a thermonuclear runaway which results in an XRB. Note that the accreted envelope above the 



4 The Eddington limit for the luminosity Lem corresponds to the maximum luminosity emitted by a stellar object 
for which the radiation pressure balances the gravitational attraction. For higher luminosities, the energy exceeding the 
Eddington limit generates a radiation-driven wind. In the case of accreting sources, the Eddington limit for the luminosity 
also puts a limit on the accretion rate MEdd- 
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neutron star is only weakly degenerate; that is, the degeneracy is lifted as the temperature increases 
at early times. For example, for a neutron star accreting fuel with Z/A=0.5, degeneracy is lifted 
(EFermi < ^thermal ~ kT) as the temperature exceeds roughly 0.3 GK. However, the extreme ejection 
velocities required for material to escape from a neutron star are never achieved by the moderate amount 
of energy released by nuclear reactions. Consequently neither appreciable mass ejection nor substantial 
envelope expansion results. 

Between individual bursts, the bulk of the observed, accretion-powered, X-ray flux arises from the 
hottest and innermost regions of the accretion disk. This persistent emission is not very sensitive to 
the nature of the compact object. Hence, the detection of thermonuclear X-ray bursts from a LMXB is 
actually the most reliable indicator of the presence of a neutron star as opposed to a black hole. 

Owing to the novel associated nucleosynthesis patterns, as well as the corresponding interest of 
experimentalists in measuring the different interactions involved, we will now focus on mixed H/He 
accretors throughout (see e.g., Refs. |24j[27] for sources exhibiting bursts in He-rich environments). 
During the thermonuclear runaway, a H/He-rich envelope is transformed to matter strongly enriched in 
heavier nuclei via the ap-process and the rp-process. These two processes involve proton- or a-particle 
induced reactions on stable and radioactive nuclei, interspersed by occasional /3-decays, depending on 
the nuclear structure of the nuclide reached by the abundance flow. Since the peak temperatures reached 
during the thermonuclear runaway may exceed 1 GK, many different nuclear interactions are likely to 
occur during the XRB. 

In this review we begin by discussing observable features of type I XRBs, highlighting characteristics 
in these phenomena that may depend upon nucleosynthesis. We then focus upon treatments that have 
been used to model XRBs and examine how the predictions of these calculations compare to one 
another and to observations. Next, we explain the roles of both theoretical and experimental nuclear 
physics in models of XRBs, and discuss nuclear physics uncertainties that either have been or remain 
to be addressed to better constrain model predictions. We conclude with suggestions to guide future 
observational, theoretical and experimental nuclear physics efforts to clarify nucleosynthesis in type I 
XRBs. 



2 Observational features of type I X-ray bursts 
2.1 What do XRBs look like? 

Most LMXBs are transient X-ray sources, indicating that the rate of accretion from the disk onto the 
compact object is not constant. LMXBs are typically discovered during periods of enhanced accretion 
rate, with persistent luminosities in the range 0.001-0.1 Lem- These periods are known as outbursts, 
and each may last for a few weeks or months. After an outburst the source returns to quiescence, 
disappearing from the X-ray sky for several months or decades. Type I X-ray bursts typically occur 
during outbursts, when the accretion rate is high enough to supply sufficient fuel for a thermonuclear 
runaway. Note, however, that some XRBs have been detected from LMXBs in apparent quiescence (6]. 

Type I X-ray bursts appear as short flashes on top of the persistent, accretion-powered emission. 
Examples of light curves are shown in Fig. [T} Depending upon the source, bursts may occur either 
quite regularly (the most regular burster being GS 1826-24) or follow a highly variable and irregular 
pattern. If the burst recurrence time is well determined, the ratio of the integrated energy released 
between X-ray bursts to the integrated energy released during the burst defines the dimensionless a 
parameter. Assuming that all accreted fuel is burned during the burst, a should then correspond to the 
ratio of energies, per unit mass, liberated via accretion and thermonuclear burning. Indeed, observed 
ratios of time- integrated persistent and burst fluxes (a 40 — 100 [IT]) match estimates based on the 
gravitational potential energy released by matter falling onto a neutron star during accretion (« 200 
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MeV per nucleon) and the nuclear energy released during a burst (~ 5 MeV per nucleon, for fuel of solar 
composition). These observations are also in agreement with detailed model calculations (see Section 



3.3.20 . 

In general, observations of X-ray bursts use detectors with medium to high time resolution. For 
example, the PCA instrument onboard RXTE could measure light curves with a time resolution of 
1 /is over an energy range of 2-60 keV, while the EPIC-pn instrument onboard the XMM-Newton 



observatory 28 has a time resolution of 30 /is over an energy range of 0.1-10 keV. These detectors also 
have medium to high spectral resolution, from E/AE m 20 — 50 for imaging cameras like EPIC, up 
to E/AE m 200 — 800 for spectrometers. This makes it possible to obtain time-resolved spectroscopy 
of X-ray bursts, which shows that during the entire burst the source spectral energy distribution is 



usually well approximated by a blackbody spectrum (but see Section 2.2). From the blackbody fit to 
the observed spectra, the colour (or blackbody) temperature Tbbody and the bolometric flux Fboi (rate of 
energy emission per unit area over the whole energy spectrum) are obtained. If one assumes isotropic 
emission from a spherical surface, the radius Rns of the emitting object (as observed at infinity) can 
then be determined using 

p _ al bbody n NS ( , s2 m 

tboi - — ~Jajji — y L + z > ( L > 

where a is the Stefan-Boltzmann constant, D is the distance from the observer to the X-ray burster, 



Mns is the mass of the neutron star and f c is the colour correction factor (see Section 2.2). The 
gravitational redshift z can be expressed as 1 + z = (1 — 2GMns/Rnsc 2 )^ 1 ^ 2 , which for a typical 
neutron star of M^s — 1-4 M Q and R^s — 10 km is (1 + z) — 1.31. Many X-ray bursters are located 
in globular clusters with well-determined distances. Measurements of radii of the emitting surface of 
X-ray bursters indicate values around 10 km |16|. 

Studies of the blackbody temperature evolution during XRBs show an initial increase in temperature 
during the increase in flux from the burst, namely, when thermonuclear burning occurs. This is followed 
by the cooling of the photosphere during the tail of the burst. This evolution occurs with a constant 
emitting area, that is, without significant expansion of the neutron star envelope (except in the particular 



cases with photospheric radius expansion, see Section 2.2). 

In general, the persistent X-ray flux F p between X-ray bursts is taken as an indicator of the accretion 
rate M, as [n] 

M — 1.33 x 10- f^V m) V. (2) 



10- 9 erg cm~ 2 s' 1 J \ 10 kpc J \IAM & ) V 1.31 / \10Jfcm 

The quantity Cboi is a correction that must be made to convert the detected X-ray flux F p to the 
bolometric flux F^i- 

Since X-ray bursts are powered by the fuel accumulated between bursts, it is expected that the 
burst rate increases with increasing persistent flux, i.e., if material piles up faster, less time is needed 
to accumulate sufficient fuel for a burst. Indeed the behaviour of some bursters is in accord with this 
expectation (see e.g., Fig. 1 in Ref. 130 1). However, in several cases, the recurrence rate instead decreases 
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as the persistent emission increases [11] . This may indicate that He is burning between bursts 31 
that the local accretion rate is actually decreasing in spite of an increase in the total accretion rate 
or that the persistent emission is in fact not a reliable indicator of the accretion rate. 

Observational characteristics of X-ray bursting sources, e.g., time between bursts, rise time of the 
light curve, burst duration, or total energy emitted, depend on various parameters such as the accretion 
rate, t he fuel composition or the neutron star spin. For accretion of solar metallicity material, Refs. 
identified several ignition regimes as a function of the adopted mass-accretion rate M: 



35 
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• M < 2 x 10~ 10 M yr _1 : Simultaneous H- and He-burning triggered by thermally unstable H- 
ignition. Here, the temperature in the accreted layer is too low for stable hydrogen burning 
between bursts. H ignites unstably in this case, triggering He burning and giving rise to a type I 
X-ray burst in a H-rich environment. 

• 2 x 10~ 10 M yr _1 < M < (4.4 - 11.1) x 1CT 10 M yr" 1 : Pure He-burning following completion 
of H-burning. Here, the temperature in the accreted envelope is higher and H burns steadily to 
He between X-ray bursts. The bursts occur as a consequence of He ignition in a pure He layer, 
resulting in fast and bright type I X-ray bursts. 

• M > (4.4 — 11.1) x 1CT 10 M yr _1 : Simultaneous H- and He-burning triggered by thermally 
unstable He-ignition. Here, accretion provides fresh H fuel at a higher rate than it is steadily 
burned. He ignites in a mixed H/He environment, resulting in type I X-ray bursts of longer 
duration. 

For sufficiently high accretion rates (M > 3.3 x 1CT 8 M yr _1 ) stable burning is expected (for fuel 
of solar composition [37]). The exact mass- accretion rate boundaries for the different regimes are still 
under investigation through detailed models. For example, for accretion of metal-deficient matter, lower 
values for the mass-accretion rates at the transition between the different ignition regimes have been 



obtained, narrowing the predicted regime for pure He-burning 34 . As well, recent simulations 37 
predict that unstable burning persists for accretion rates higher than those inferred from observational 
constraints, calling into question other assumptions in the models (e.g., the composition of the accreted 
material). 

Since one source may evolve through different accretion rates during an outburst, in principle the 
different burning regimes could be observed in the same source. In several bursters, the transitions 
between burning regimes have nicely shown up through the observation of changes in burst durations 



30 ,38 39 , but exceptions are abundant. Many sources show the exact opposite behaviour from what 



is expected: burst durations decrease as the apparent accretion rate increases for sources thought to be 



passing from the second to the third regime above 40 -44 



2.2 Photospheric radius expansion (PRE) bursts 

In some bright type I X-ray bursts, the luminosity reaches the Eddington limit and causes the expansion 



of the envelope. This is illustrated in Fig. |2| for a type I burst observed from the Rapid Burster 45 
A first indicator of the expansion of the envelope is a flat-topped light-curve, with a nearly constant 
bolometric luminosity during the envelope expansion (top panel). Confirmation of the expansion is 
provided by time-resolved spectroscopy of the constant bolometric luminosity interval, showing the 
increase (expansion) and decrease (contraction) of the photospheric radius (middle panel) in parallel 
with the decrease and increase of the effective temperature (bottom panel). This is followed by the final 
cooling of the neutron star envelope. See also, e.g., Fig. 3.4 from Ref. |46| for examples of both PRE 
and non-PRE bursts from a single source. 

In some extreme cases, the temperature drop is large enough to shift the X-ray emission to energies 
below the range of sensitivity for typical instruments, causing an apparent double-peak in the detected 
X-ray flux. This effect should not be confused with intrinsic double or multiple peaks of the bolometric 



flux (see Section 2.4) 



Observations of these PRE bursts probe the still-unknown equation of state (EOS) of the neutron 
star through the determination of the neutron star's mass and radius. Several models for the EOS have 



been developed, with and without condensates, and some including strange quark matter 47-49 . As 
each EOS predicts a particular relation between the mass and radius of the neutron star, these EOS 
models can be tested using observational constraints on masses and radii. For a PRE burst located at 
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a known distance, the neutron star mass can be determined from the observed bolometric luminosity if 
one assumes that the maximum observed luminosity corresponds to the Eddington limit. The neutron 
star radius can then be inferred from spectral fits during the cooling tail of the burst if one assumes the 



emitting area to be the entire neutron star surface 49 52 



A number of details play a role in the determination of the neutron star radius and mass from PRE 
bursts, and the composition of the envelope is a critical element in these calculations. The hydrogen 
mass fraction X is needed to determine the electron scattering opacity K e — 0.2(1 + X) cm 2 g _1 , which 
in turn is used to determine the Eddington limit for the luminosity (and hence the neutron star mass). 
In addition, the composition of the atmosphere determines the colour correction factor: the neutron 
star spectrum is known to differ from a perfect blackbody. This is because electron scattering in 
the hot neutron star atmosphere suppresses part of the emission 53 . The colour correction factor 
fc = Tbb 0< iy/T e ff relates the effective temperature T e ff to the colour or blackbody temperature T^ d y 
obtained from the spectral fit. The quantity f c affects the determination of the emitting area, and hence 
the neutron star radius, since the emitting area is obtained from the observed flux and temperature. 
Detailed calculations using atmosphere models 54 56 have determined colour correction factors for 
metallicities up to solar; the smallest colour correction factors were obtained assuming solar abundances 
because the absorption edges of Fe move the neutron star spectrum closer to a blackbody spectrum. 
However, in Eddington-limited bursts, the photosphere of the neutron star is expected to contain high- 
metallicity ashes from the thermonuclear burning of H and He deep in the envelope 57 58 . This is in 



agreement with results from hydrodynamic simulations, which show that convective mixing approaches 
the photosphere when the flux reaches the Eddington limit 59 61 . Increased atmospheric metallicity 



would imply more absorption features in the final emerging spectrum and even smaller colour correction 
factors. Atmosphere model calculations adopting metallicities larger than solar are therefore urgently 
required to examine the impact of smaller colour correction factors on neutron star mass and radius 
determinations using PRE bursts. 

In addition to the uncertainties introduced by the differences between an ideal blackbody spectrum 
and the actual emission from the neutron star atmosphere during a XRB, the final observed spectrum 
may be altered by the interaction of the emerging spectrum with the accretion flow or the hot corona 
surrounding the neutron star (e.g., through reflection, partial occultation, Comptonization) . The re- 
sulting uncertainties seem to be sufficiently small 51 , 52 to still allow one to determine neutron star 
masses and radii with sufficient precision to distinguish between different neutron star EOS models. 

There are therefore two ways in which the ashes of thermonuclear burning in XRBs may influence the 
determination of neutron star masses and radii from the observation of PRE bursts: via the opacities 
of the expanding layer (affecting the value of the assumed Eddington limit L Edd ) and via the inter- 
pretation of the observed spectrum (through the deviations of the actual atmosphere spectrum from 
a perfect blackbody). Note that other quantities characteristic of neutron stars may also be affected 
by the composition of burst ashes. These include thermal 62,63 , radiative 64 , electrical 63 , 65 and 



mechanical 66 67 properties of the neutron star crust, which are important for understanding, e.g. 
the evolution of the neutron star magnetic field or the persistent emission between bursts. 



2.3 Burst oscillations: burning in confined regions? 



Oscillations in the light curves of X-ray bursts, with frequencies between about 250 and 600 Hz 46 



were first reported by Strohmayer et al. (1996) p38] and have now been identified in about 25% of all 
bursting sources. In bursts from LMXBs, oscillations during the rise time may be consistent with the 



spreading of a hot spot on a rotating neutron star 69 . Oscillations during the burst tails, however, 
exhibit a drift in frequency, increasing by a few Hz and approaching an asymptotic value as the burst 
progresses. The asymptotic frequency of a burst seems to be a characteristic attribute of an individual 



source, being stable from burst to burst 70 . Oscillations have also been observed in bursts from 
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accreting millisecond pulsars 71 



no noticeable drift during the burst tails 72 , 73 



at frequencies corresponding to the known spin frequencies and with 
This suggests that a confined radiating region would be 



the most straightforward way to account for the oscillations, although more sophisticated explanations 
are clearly needed to account for the frequency drifts in bursts from LMXBs. Indeed, altogether different 
underlying mechanisms may be responsible for the burst oscillations in the two scenarios. Suggested 
mechanisms include surface wave modes in the neutron star ocean and atmosphere 74 77 , Coriolis 
force containment [78] , and magnetic confinement of the fuel 79 . While the magnetic field in most 
LMXBs is likely too small (~10 8 G) to confine the fuel, the magnetic field could be enhanced during 
bursts 80 



2.4 Multi-peaked X-ray bursts 



Several bursters have exhibited multi-peaked X-ray bursts without evidence of photospheric expansion, 
where both the bolometric luminosity and the temperature show multi-peaked profiles. These include 
double-peaked bursts from 4U 1608-52 [IT], GX 17+2 [82] and 4U 1709-267 



83 , and about a dozen 



double-peaked and two unique triple-peaked bursts from the well-studied source 4U 1636-53 11,41,84 



Examples of double-peaked light curves are shown in Fig. [TJ The separation of the two peaks in 
the bolometric light curve is typically ^3-5 seconds, with the flux decreasing to values beween 20% 
and 50% of the maximum peak flux during this interval (TTJ [86j [87] . Time-resolved spectral analysis 
shows that the emitting area increases steadily from the start of the bursts until the beginning of the 
cooling tail 86,87 . As well, a study of double- and triple-peaked bursts from 4U 1636-53 showed that 



all multi-peaked bursts occurred during periods of high accretion rate 39 



Bhattacharyya & Strohmayer 86,87 suggested that aspects of the observed double-peak features 



could be explained through the propagation of a burning front ignited either at the neutron star rota- 
tional pole or at high latitudes. The propagation towards the equator would be stalled at low latitudes 
as it would advance against the flow of matter transferred onto the equatorial plane from the accretion 
disk. This could explain the cooling after the first peak. After crossing the equator, the burning front 
would advance rapidly again towards the opposite pole, explaining the second peak (but see Ref. [89]). 
In addition, ignition at the rotational pole would imply the absence of oscillations in the observed 
emission, as measured in some cases [86] . Ignition at high latitudes (but not at the pole) could explain 
the millisecond oscillations detected in the first peak for 4U 1636-53, while the spread of the burning 
flame throughout the entire surface would explain the disappearance of the oscillations during the sec- 
ond peak [87] . Other explanations include two-step energy generation due to shear instabilities in the 
fuel |90| or scattering of the X-ray emission by material evaporated from the disk 91 , but these models 
do not reproduce the observed double-peaked burst temperature and radius profiles 86 . A nuclear 



waiting point impedance in the thermonuclear reaction flow (see Section 3.3.2) was suggested by Fisker 
et al. [92], but this has difficulties in explaining the large dips observed between the two peaks 86 . 



2.5 Absorption features in X-ray bursts 

The detection of absorption lines in the neutron star atmosphere is potentially the most direct and 
powerful tool for probing both the neutron star equation of state (via determination of the gravitational 
redshift z), and the products of nucleosynthesis during XRBs (via the determination of abundances in 
the neutron star atmosphere, especially for PRE bursts |57|). The latter consideration is particularly 
relevant given that no material has been observed to be ejected during XRBs. The main observational 
difficulty with the detection of absorption features in XRBs is the short duration of the bursts: long 
exposure times are required to accumulate data with enough signal to noise. Even for the most regular 
bursters such as GS 1826-24, KS 1731-260 and IM 0836-425, with bursts of 20-40 s in duration occurring 
every 2-3 hours, only ^0.3% of the total observation time is spent measuring the emission from the 
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bursts. As well, line broadening effects can exacerbate the difficulties in identifying and interpreting 
93H95] (but see also Ref. pi). 



absorption features 
Cottam et al. 



96 



reported the first detection of gravitationally-redshifted absorption lines, in the X- 
ray spectrum of EXO 0748-676. These results were serendipitously obtained as EXO 0748-676 had been 
chosen as a calibration source for the XMM-Newton observatory. As a result, it was observed for 335 ks 
with the onboard Reflection Grating Spectrometer; such long periods for spectroscopic observations are 
not normally possible with X-ray observatories. A total of 28 type I bursts occurred during this time, 
with an effective exposure time to the X-ray burst spectrum of only 3200 s. This allowed the authors 
to identify several absorption features, including lines from Fe XXV, Fe XXVI and O VIII that had 
been gravitationally redshifted by z = 0.35. However, these measurements were not confirmed during 

In addition, detailed 
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a second, longer observation of the same source with the same instrument 
calculations of the expected emission spectrum from the neutron star atmosphere of EXO 0748-676 |98] 
have cast serious doubts on the ability to identify absorption lines. A long observation of the regular 
bursting source GS 1826-24 also failed to detect any absorption features 
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Exciting results have recently been reported using RXTE that indicate the observation, albeit with 
limited spectral resolution, of absorption edges from Fe-peak elements in photospheric radius "super- 
bursts 



expansion 



58 . Such efforts are highly encouraging and future studies using instruments with 



higher spectral resolution (such as XMM-Newton or Chandra |100| ) should be a priority. 



3 Modeling type I X-ray bursts 



3.1 Defining the mechanism 

Within a year of the serendipitous discovery of XRBs [l][5] more than 20 bursting sources were discovered. 
This included the enigmatic Rapid Burster 14 , 101 , clearly at odds with the general behaviour exhibited 



by all the other bursting sources (see Section lj). Any model attempting to characterize XRBs had to 
account simultaneously for the intermediate to long recurrence periods inferred from the majority of 
the sources as well as for the extremely short recurrence periods reported for the Rapid Burster. 

A first series of models, in which bursts were driven by instabilities in the accretion flow onto 
a compact, stellar corpse (a white dwarf, a neutron star, or a black hole) was proposed by several 
groups. Clearly, a compact object was required to host the explosion to guarantee that large amounts 
of gravitational energy were released as X-rays by matter falling into its deep gravitational well. The 
exact nature of the underlying compact object was a matter of debate, and in fact one of the early 
models assumed accretion onto a supermassive black hole (> 100 M 102] ) . A breakthrough was 
achieved through observations of the bursting source 4U 1724-30 with the OS 0-8 satellite 103 : the 
spectral evolution of one of its long bursts, best fitted with a blackbody spectrum with kT « 0.87 - 
2.3 keV, suggested a much smaller source (i.e., either a neutron star or a stellar black hole). Assuming 
a source distance of ~ 10 kpc, Refs. 104, 105 inferred a blackbody radius for the source of ~ 10 km. 



Other features such as spectral softening during the decay of type I XRBs, the harder X-ray spectra of 
bursting sources as compared with most of the X-ray transients hosting black hole candidates, as well 
as the masses inferred from those same systems, pointed towards neutron star primaries 106 . These 



pioneering simulations considered effects such as instabilities in the interaction between the accreted 



matter and the neutron star magnetosphere 107 114 , flare-like eruptions of magnetic energy generated 
and released through a toroidal field in the accretion disk 115 , convective-driven instabilities in the 



disk [116 . 117 , or thermal instabilities driven by Compton heating of the accretion flow 118 



A second series of calculations relied on thermonuclear explosions in the layers accreted onto a 

It was found that the characteristic timescale of a thermonuclear runaway 



neutron star 



1243 120 122 



on top of a neutron star was < 1 s [Tl9] , in rough agreement with the timescale of the X-ray emission 
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reported from some bursting sources. It was clear, however, that the quick succession of flashes observed 
from the Rapid Burster could not be accounted for by this thermonuclear runaway model since not 
enough fuel could be piled up on the surface of the neutron star in such short time intervals. Progress 
was made through the subsequent observation of two different types of bursts associated with the Rapid 
Burster, with different recurrence periods 15 123 . The classification of type I and type II bursts (the 



former associated with thermonuclear flashes, the latter to accretion instabilities) was then suggested. 



3.2 The thermonuclear runaway model 



Rosenbluth et al. 124 were the first to quantify the nuclear energy released from the fusion of H-rich 



material on accreting neutron stars. It was later concluded 119 125 that under certain conditions such 



nuclear fusion episodes may trigger thermonuclear runaways. The connection between unstable nuclear 
burning on the surfaces of neutron stars and typ 
(for He-driven bursts) and Maraschi and Cavaliere 



of increasingly detailed numerical simulations 35 , 64 



predictions and observations were made by Refs. 120 
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burning on the surfaces of neutron stars and type I XRBs, as envisaged by Woosley and Taam 12 

(for H-driven bursts)^] paved the way for a suite 



127f 136 . Initial efforts to compare theoretical 



121 



The gross observational features of type I XRBs were qualitatively reproduced by early studies, 
already suggesting the possibility of constraining neutron star properties (e.g., masses, radii, equation 
of state) through models. Early dimensional analysis work 120 121 revealed that for bursts driven 



by He-burning, one could expect recurrence periods of about 10 000 s, accreted envelope masses of the 
order of ~ 10 21 g, and an overall energy release of < 10 39 erg per burst. Later predictions included peak 
luminosities of ~ 10 38 erg s -1 (with peak blackbody temperatures of 3 x 10 7 K), light curve rise times 
of ~ 0.1 s, and light curve decay times of ~ 10 s, through consideration of the transport of nuclear 
energy from deep in the envelope to the outermost neutron star layers. Another important observational 



constraint matched by thermonuclear models of type I XRBs was the a parameter (see Section 3.3.2), 
the ratio between time-integrated persistent and burst fluxes. 

The first detailed numerical models of XRBs can be found in Joss (1978) |127 for several sets 



of mass-accretion rates and neutron star central temperatures. (The central temperature determines 
the initial surface luminosity assumed in a model.) These parameters, together with the neutron star 
mass, equation of state and the metallicity of the accreted material, represent key ingredients in XRB 
modeling. Joss's work quantitatively confirmed the values derived from dimensional analysis. As well, 
it stressed that a burst consumes virtually all the accreted fuel (likely synthesizing Fe-peak elements) 
and that the energy released during the explosion is essentially emitted as X-rays from the neutron star 
photosphere. Further numerical simulations 128 revealed that for highly magnetized neutron stars 



(B > 10 12 G), matter transferred from the stellar companion would be funneled onto the neutron star 
magnetic poles, enhancing the local accretion rates in those regions by a factor of ~ 10 3 . Ignition 



regimes as a function of mass-accretion rate were later defined 35 , 36 for the accretion of material of 



solar metallicity (see Section 2.1). 



General relativistic corrections were first incorporated into XRB models by Refs. 129| 130[ 132 
more consistently by Ref. 



135 



and 



These are important for the interpretation of observational properties 
of X-ray bursts. For example, previously predicted peak luminosities and burst recurrence times were 
reduced and increased by a factor of 1 + z, respectively, where z is the gravitational redshift (see Section 



2.1). 



5 According to Ref. 16 , the earliest suggestion of the thermonuclear origin of type I XRBs was made by L. Maraschi 
while visiting MIT in February 1976, at the time when the central region of the Milky Way was being observed with the 
SAS-3 satellite fl26l. 



10 



3.3 Model predictions of nucleosynthesis in XRBs 



3.3.1 Overview 



Before discussing results from XRB models, one must carefully distinguish between the different nu- 
merical approaches used in the calculations. One category of models includes parameterized one-zone 
calculations 



35 63 137 141,153 



These simple prescriptions relate the history of the neutron star's 
accreted envelope with the time evolution of the temperature T and density p in a single layer of the 
envelope. Such thermodynamic quantities are often determined through semi-analytical models, or 
occasionally correspond to T-p profiles directly extracted from 1-D hydrodynamic models. This ap- 
proach, while representing an extreme oversimplification of the physical conditions governing neutron 
star envelopes, was used extensively to overcome the time limitations that arose in computationally 
intensive hydrodynamic calculations. For similar reasons, one-zone calculations have been used recently 



to estimate the impact of nuclear uncertainties on the final XRB yields 142 145] (see Section 4.3). 



At the moment, state-of-the-art XRB nucleosynthesis calculations rely on 1-D hydrodynamic models 



(see Refs. 59 61 ,136 146 148 , and references therein). The underlying assumption in all these models 
is spherical symmetry. This simplifying hypothesis demands that the explosion occur uniformly over a 
spherical shell. Moreover, spherically symmetric models are limited by their treatment of the manner 
in which thermonuclear runaways are initiated (presumably as point-like ignitions) and propagate. 
This particularly affects convective mixing, a key ingredient in deflagrative, explosive scenarios such 
as type I X-ray bursts that can only be accurately modeled in three dimensions. Multidimensional 
XRB models are computationally challenging and only preliminary simulations of specific aspects of 
the explosions (such as flame propagation 78 or the early convective stages preceding ignition 230 1) 
have been conducted to date. Hence, no reliable nucleosynthetic predictions have been reported using 
multidimensional simulations. 

The relevant nuclear reaction sequences accompanying XRBs were initially identified in the 1980s- 
90s 63 137 150 151 . Because of computational constraints, early attempts to predict the nucleosynthe- 



sis expected during XRBs were performed using limited nuclear reaction networks, truncated near Ni 
(and using a 19-isotope network) [136 147 148] ; Cu (181-isotope network) |139| ; Kr (274-isotope |138| 
or 463-isotope networks |140| ); Cd (16-isotope network) [152]; or Y (250-isotope network) |137 . Efforts 
to link state-of-the-art hydrodynamic models with both accurate input physics and large nuclear reac- 
tion networks were clearly needed but required advanced generations of computers not available until 
recently. 



3.3.2 Results from models with large nuclear reaction networks 

With a neutron star as the underlying compact object hosting a thermonuclear explosion, temperatures 
and densities in the accreted envelope reach T > 10 9 K and p ~ 10 6 g cm -3 . Detailed nucleosynthesis 
studies under such conditions require the use of networks with hundreds of isotopes and thousands 



of nuclear interactions. Pioneering efforts in this regard were made by Schatz et al. 63 153] , who 
carried out XRB nucleosynthesis calculations with a network containing more than 600 isotopes (up 
to Xe). Unfortunately, those calculations relied on a one- zone approach and so inadequate attention 



was given to the stellar model (and to convection in particular). Koike et al. 141 also performed 
detailed one-zone nucleosynthesis calculations, using temperature and density profiles extracted from 
spherically symmetric hydrodynamic models. These profiles were linked to a 1270-isotope network 
extending up to 198 Bi. Despite the limitations in the adopted approach, these works confirmed that 
the main nuclear reaction flow is driven by the rp-process (rapid proton-captures and /3 + -decays), the 
3a-reaction, and the ap-process (a sequence of (a,p) and (p/y) reactions). They revealed that the 
abundance flow proceeds far away from the valley of stability, merging with the proton drip-line beyond 
A = 38. They also stressed the need to rely on huge networks, since the main nuclear path extends 
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up to the SnSbTe mass region |153| or even beyond (e.g., the nuclear activity in Ref. 141 reached 
126 Xe). It is clear that the location of the nucleosynthetic endpoint in XRB simulations depends upon 
the astrophysical parameters adopted. The specific location of the maximum possible endpoint in the 
most extreme bursts is still a matter of debate, since recent experimental studies 154 have revealed 
difficulties in reaching the SnSbTe-mass region. 

Only very recently has it been possible to couple hydrodynamic models and detailed nuclear reaction 
networks. Efforts include Rcfs. [37, 60, 92, 155, 156], with networks of ~300 isotopes (up to Te); Rcfs. 
59,157 , with networks of 1300 isotopes and an adaptive network; and Refs. 61 , 158| , with networks 
of 325 isotopes (up to Te). The first set of calculations (37][60|[92||155 156 relied upon the 1-D, 
fully relativistic hydrodynamic code AGILE, with 129 numerical shells and convective transport based 
on mixing-length theory. In turn, the second set 59, 157] was performed with the 1-D, Newtonian 



hydrodynamic KEPLER code, with a variable number of shells (up to 1000 for some models), time- 
dependent convective transport based on mixing length theory with overshooting, and thermohaline 
mixing. Finally, the third set of calculations 61 , 158 used the 1-D, Newtonian hydrodynamic code 



SHIVA, with up to 200 numerical shells and convective transport based on time-dependent mixing-length 
theory. Overall, a number of XRB observables have been reproduced by current hydrodynamic models. 
This includes recurrence periods, peak luminosities and a values (see, e.g., Heger et al. (2007) [33] for a 
comparison between model predictions and the regular bursts from GS 1826-24). Nonetheless, certain 
aspects associated with the observed shapes of the light curves, such as the rise times, the presence of 
bolometric light curves with multiple peaks, or the incredible diversity in light curve shapes observed 
from individual sources, have not been fully understood yet. 

Here we summarize results obtained by the Barcelona group (6l] using models of 1.4 M neutron 
stars accreting matter of different metallicities at a rate of 1.75 x 10~ 9 M yr -1 (0.08 Mem)- The initial 
luminosity of the neutron star, at the onset of accretion, is L initia i = 4.14 L Q . We outline the most 
striking features of the main nuclear activity during X-ray bursts and compare these results with those 
from simulations performed by the Santa Cruz 59 and Basel 60 groups. Interested readers can find 



more detailed accounts of the dominant flows at different stages of the bursts in Refs. 59 -61 



For models accreting material of solar metallicity, results include burst recurrence times of r rec ps 5 
hr, peak temperatures T pea k of ~ 1.1 — 1.3 GK, peak luminosities L pea k = (4 — 8) x 10 38 erg s -1 and 
ratios between persistent and burst luminosities in the range of a = 30 — 40. While these values are 
qualitatively in agreement with those inferred from the bursting sources GS 1826-24 (r rec = 5.74 ± 0.13 
hr; a = 41.7 ± 1.6), 4U 1323-62 (r rec = 5.3 hr; a = 38 ± 4), or 4U 1608-52 (r rec = 4.14 - 7.5 hr; 
a = 41 — 54), the simulations systematically yield larger recurrence times and peak luminosities (and 
hence, lower a) than those reported by the Santa Cruz group 59 for their model ZM. The twelve 
bursts computed in Ref. 59 are characterized by recurrence times of about ps 2.7 hr, peak luminosities 
of L peak pa (1.5 — 2) x 10 38 erg s -1 , and ratios between persistent and burst luminosities of a ~ 60 — 65. 
We note that models from the Basel group 



60 



give lower L veak and somewhat larger a values than 
models from both the Barcelona and Santa Cruz groups. 

The time evolution of temperature and density at the ignition shell during the third bursting episode 
of the model from the Barcelona group is shown in Fig. |3j top and middle panels. Note the modest 
expansion of the envelope during a burst (Fig. |3j lower panel), which takes place at nearly constant 
pressure. The moderate T veak achieved in this model restricts most of the nuclear activity to A ps 60 
(mainly 60 Ni and 64 Zn |61|) and no large concentrations of material in the SnSbTe-mass region were 
found, in agreement with the results obtained by the Santa Cruz group |59|. The mean composition of 
the envelope at the end of the fourth burst in this model is illustrated in the top panel of Fig. [4] (filled 
circles); the corresponding light curve is given in the lower panel of Fig. [4] (black line). This model 
also yields a very small post-burst abundance of 12 C, below the threshold amount required to power 
superbursts (X( 12 C) min > 0.1, see Ref. |20|). 



Results from a second set of models, identical to those described above but considering accretion 
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of metal-defficient matter (Z = Z /2O = 0.001), include longer recurrence times of r rec w 9 hr, T pea k 
of about 1.3 — 1.4 GK, L peak « 10 38 erg s _1 and ratios between persistent and burst luminosities of 
a ~ 20 — 30. These values are comparable to those measured in the XRB sources 1A 1905+00 (r rec = 
8.9 hr), 4U 1254-69 (r rec = 9.2 hr), or XTE J1710-281 (r rec = 8.9 hr, a = 22-100). The longer exposure 
to high temperatures (driven by a slower decline phase) causes an extension of the main nuclear path 
towards the SnSbTe-mass region, with final abundances at the end of the fifth burst dominated by the 
presence of 105 Ag, 64 Zn, 104 Pd and 68 Ge. The mean composition of the envelope at the end of the fifth 
burst in this model is illustrated in the top panel of Fig. [4] (open squares); the corresponding light curve 
is given in the lower panel of Fig. [4] (grey line). Similar to the solar metallicity model, little 12 C remains 
at the end of a burst. 

The two models from the Barcelona group discussed above differ only in the metallicity of the 
accreted material. The nuclear activity in these two models is compared in detail in Figs. [5j [6j [7]and|8j 
These plots show the most abundant species (X > 10 -4 ) and dominant reaction fluxes at T pea k, within 
the ignition shell of the fourth computed burst in each model. Figs. [5] and [6] reveal the existence of 
waiting point nuclei: unstable nuclei (e.g., 30 S, 56 Ni, 60 Zn, 64 Ge, 68 Se, 72 Kr, 76 Sr, 80 Zr, 89 Ru, etc.) where 
the abundance flow is temporarily delayed, allowing material to accumulate. The manner by which each 
waiting point nuclide is eventually depleted depends on the nuclear structure in the immediate mass 
region of this species: some are destroyed by /3-decay, some by (a,p) reactions, and others by sequential 
(rather than direct) two-proton capture (see Section 4.1). The bursts in both models are eventually 
quenched by the exhaustion of H and He fuel rather than through any significant envelope expansion (see 
Fig. [3j lower panel). We note that a typical range of metallicity for the accreted material is difficult to 
determine. Whereas solar metallicities may be regarded as representative (e.g., a thorough comparison 
between observations of the X-ray burster GS 1826-24 and models with different metallicities favoured 
the accretion of material of solar metallicity [33]), the identification of many bursting sources in globular 
clusters implies that lower metallicities should also be expected. 

Comparison of results from the solar and low-metallicity models suggests that the smaller the metal 
content of the accreted material, the larger the recurrence time (and the smaller the a). In addition, 
explosions in metal-deficient envelopes are characterized by lower peak luminosities and larger light 
curve decline times. Similar trends are reported by the Santa Cruz group 59 when comparing models 
accreting solar and low metallicity material. This is in spite of any systematic differences between the 



models, as mentioned above. The most striking difference between the models of Ref. 59 and Ref. 61 
concerns the larger effect played by the metallicity of the accreted material in the models computed 
by the Barcelona group [6l]. Indeed, Ref. [59] explained the moderate effect they found as due to 
"compositional inertia" washing out the influence of the initial metallicity. This effect refers to how 
properties of bursts may be sensitive to the fact that accretion occurs onto the ashes of previous bursts, 
minimizing the impact of the actual metallicity of the accreted material after sufficient bursts have been 



calculated. Similar conclusions were reached by Ref. 60 . The origin of these discrepancies (as well 
as other discrepancies in peak luminosities, recurrence rates or a values for a common metallicity) is 
not yet clear. Although XRB properties depend only weakly upon the neutron star mass (or surface 
gravity), some of the differences outlined between the three studies may be attributed to the combined 
effects of different assumptions for the adopted neutron star size and to differences in the input physics 
(e.g., nuclear reaction network, opacities, treatment of convective transport). 

Aside from any metallicity-dependent effects, it is evident from Fig. [4] that different nucleosynthesis 
patterns are associated with different characteristic burst light curves. This is seen even within sequences 



of bursts calculated at a given metallicity (see, e.g., Figs. 19 and 20 in Ref. 61 , or Figs. 13 and 15 



in Ref. |59|). This relationship is of vital importance since the products of nucleosynthesis during 
an XRB have not been directly observed in a straightforward, reliable manner (see Section 2.5). As 
such, from the point of view of nuclear astrophysics, experimental nuclear physicists cannot rely solely 
on exploring the effect of a nuclear laboratory measurement on XRB nucleosynthesis; one must also 
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examine the effect of the measurement on the principal observable of an XRB, namely, the burst light 
curve. This is in stark contrast to the cases of classical nova explosions or supernovae, in which the 
products of nucleosynthesis are ejected into space and can be observed through spectroscopy or other 
means. Fortunately, the link between nucleosynthesis in XRBs and the predicted light curves provides 
an essential guide to help judge which nuclear physics uncertainties may most significantly affect XRB 



light curves (see Section 4.3). 

The effect of the adopted mass-accretion rate on burst properties has recently been investigated 
by Ref. |59j. Models of 1.4 M neutron stars (L iniUa i = 4.14 L ) accreting matter at 3.5 x 1CT 10 M 



yr _1 (0.02 M^dd) and at 1.75 x 10~ 9 M yr _1 (as mentioned above) were computed for two different 
metallicities (solar and 5% solar). For the models accreting material of solar metallicity, those with 
M = 1.75 x 10~ 9 M exhibited combined H/He flashes, while bursts from models with M = 3.5 x 10~ 10 
M yr -1 experienced ignition in a pure, H-free, He layer, yielding briefer, brighter light curves with 
shorter tails, very rapid rise times (< 0.1 s), and ashes corresponding to elements below the Fe region. 



This is in agreement with the discussion in Section 2.1 



With regard to bursts featuring bolometric luminosities with double-peaked profiles, Fisker et al. 92 



report that some of the observed features (see Section 2.4) could be explained, in part, through a waiting 
point impedance in the nucleosynthesis accompanying an XRB. For suitable accretion rates, the first 
peak could be caused by the rapid ignition of He accumulated at the base of the envelope, following 
the burning of H through the CNO cycle. Energy released from this He burning region would stimulate 
thermonuclear ignition in other shells where H may not have been fully converted to He. The abundance 
flows in these higher-lying shells may then determine whether or not a second peak is produced. For 
example, if the flow encounters a waiting point nucleus, the nuclear energy generation rate may decrease, 
possibly causing the surface luminosity to also decrease. Once the flow proceeds again, beyond the 
waiting point, the surface luminosity may increase again. Fisker et al. suggest that waiting points 
at 22 Mg, 26 Si, 30 S, and 34 Ar may be of consequence in this context, and that reaction rates to better 
determine how these nuclei are eventually destroyed (e.g., through (a,p) reactions) should therefore be 
measured. 

Finally, a long-debated aspect of the nucleosynthesis accompanying XRBs is the potential role of 
these phenomena as nuclear factories of the elusive light p-nuclei 92,94 Mo and 96,98 Ru 



57, 153, 159 



161| (for a review on the production of p-nuclei, see Arnould and Goriely 162| ). Matter accreted 
onto a neutron star of mass M NS ~ 1-4 M and radius R^s ~ 10 km releases GM NS fn p /R NS ~ 
200 MeV per nucleon, whereas only a few MeV per nucleon are released from thermonuclear fusion. 
Ejection from a neutron star is therefore energetically unlikely, as confirmed by all recent hydrodynamic 
simulations 59 6l]p Even though radiation-driven winds during photospheric radius expansion (see 
Section 2.2) may lead to the ejection of a tiny fraction (i.e., ~1%) of the envelope 57, 164 167 , it has 



not yet been rigorously analyzed whether such winds may contain material synthesized during the burst. 
Furthermore, in one-zone nucleosynthesis calculations the chemical species are assumed to represent the 
whole (by construction, chemically homogeneous) envelope. In multi-zone hydrodynamic simulations, 
however, the abundances of many species, including these p-nuclei, decrease dramatically towards the 
outer envelope layers because of limited convective transport. Unfortunately it is these layers that 
are most likely to be ejected by radiation-driven winds. The predicted overproduction factors in these 
outer regions [61] are several orders of magnitude smaller than those required to account for the origin of 
Galactic light p-nuclei [161 , in sharp contrast with the results reported from one-zone calculations [153] . 
We conclude that XRBs are unlikely to be dominant contributors to the Galactic abundances of p-nuclei, 



6 The fact that no mass escapes the computational domain, even in Lagrangian frameworks (i.e., simulations in which 
the computational grid is attached to the fluid), allows one to simulate sequences of repeated bursts (see pioneering 



models by Refs. 136 163] ). This has opened up new possibilities to study the long-term evolution of bursting sources 



and to identify asymptotic regimes for suitable combinations of mass-accretion rates and metallicities (see discussion in 
Ref. [59]). 
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according to current models. 



4 Nuclear physics in type I X-ray bursts 
4.1 Reaction rate formalism 



We have already discussed in Section |3.3.2| some general features of type I X-ray burst nucleosynthesis. 
The interplay of charged-particle reactions and /3-decays at temperatures and densities characteristic of 
XRBs crucially defines the overall abundance flow and nucleosynthesis. Instead of a smoothly continuous 
abundance flow from lighter to heavier nuclear species we should envision a flow in very fast spurts (<1 
s), each delayed significantly at the "waiting point nuclei" where the capture of successive charged- 
particles is inhibited by strong reverse photodisintegration reactions. From these general considerations 
we can deduce the nuclear quantities generally needed for studying the nucleosynthesis in type I XRBs. 

The basic quantity required is the thermonuclear reaction rate per particle pair, (crv), defined as 
the number of nuclear reactions taking place in the stellar plasma per unit volume and per unit time. 
Specifically, the reaction rate represents the integral over all kinetic energies E, of the product of the 
Maxwell-Boltzmann distribution and the nuclear reaction cross section a . It is given by 

where m is the reduced mass, is Avogadro's number, k is the Boltzmann constant, and T is the 
stellar temperature. There are a number of strategies for estimating the thermonuclear reaction rate, 
depending on the nuclear reaction mechanism. When the cross section varies smoothly with energy 
and is devoid of any structure, the reaction most likely proceeds via a very fast (i.e., less than 10~ 20 
seconds), single step mechanism. On the other hand, if the incoming projectile is absorbed by the target 
nucleus, sharing all its energy with the target nucleons over a much longer period of time (i.e., about 
1CT 16 seconds), the cross section most likely will exhibit a sharp peak, called a resonance. In reality 
many resonances with different energies and strengths (i.e., the integrated resonance cross section) will 
contribute to the reaction rate. In the most favourable case, the total reaction rate can be derived 
from experimental nuclear data. Significant advances have been made recently for evaluating such 
"experimental" reaction rates using a Monte Carlo procedure in order to estimate statistically rigorous 



values 168 , 169 . On the other hand, when no nuclear data are available, the reaction rate must be 
estimated with the aid of nuclear reaction models, for example, the statistical (Hauser-Feshbach) model 
(see below), or the nuclear shell model (see e.g., 170| ). In either case, reaction rates derived under the 



assumption that all target nuclei are in their ground states ("laboratory rates") must, in general, be 
corrected for thermal excitations of the target nuclei in order to estimate "stellar rates" . The required 
correction, or "stellar enhancement factor" (SEF), can also be derived from the Hauser-Feshbach model 
(see, for example, Refs. 171 , 172] ) . For the specific case of type I X-ray bursts, the peak temperatures 



achieved (< 2 GK) result in small predicted correction factors for most nuclear reaction^] 

It is often not appreciated that under certain special conditions the thermonuclear rate of a given 
reaction is irrelevant, despite the fact that the main abundance flow passes through this nuclear link. 
This situation may arise, for example, at sufficiently high temperatures when the abundance flow 
encounters a waiting point with a very small (less than a few hundred keV), or even negative, proton 
separation energy (i.e., the minimum energy required to initiate proton emission). In this case, an 
equilibrium between the rates of the forward reaction and corresponding reverse photodisintegration 
is quickly established. Clearly, once equilibrium has been established for a given pair of nuclei A and 

7 For a recent discussion of "stellar enhancement factors" , and the related but distinct "ground state fraction of the 



stellar rate" , see Rauschcr et al. 184 . 
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B, the abundance flow and nucleosynthesis near the waiting point become independent of the reaction 
rates linking these two species. To be more specific, consider nuclides A, B and C located near the 
proton dripline, where an equilibrium A + p -H- I? + 7 is established. Furthermore, assume that species 
B may only be depleted via another reaction, B + p — > C + 7. The nucleosynthesis depends then on 
the relative magnitude of the decay constants for the competing /5-decay and the sequential two-proton 
capture of the waiting point nucleus A. The latter decay constant is given by 



168 
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(4) 



where N p is the proton number density, h is Planck's constant, is the reduced mass of A + p, the 
ji are the nuclear spins, and the G 1 } 01 "" 1 are the normalized partition functions (i.e., a factor expressing 
the thermal distribution of nuclear levels in a given nuclide). The last two factors in Eq. (3) contain the 
most important quantities: the decay constant for sequential two-proton capture depends on the Q- value 
(i.e., the nuclear energy released), but not the reaction rate, of the A(p,j)B reaction and, in addition, 
on the reaction rate of the subsequent B(p, r y)C reaction. Because the Q-value enters exponentially in 
the above expression it must be known with good precision, say, better than a few tens of keV 
Particularly important are the Q-values for the (p,7) reactions on the waiting point nuclei 64 Ge 
and 72 Kr 



145 



68 



Sc 



168 . Many of the Q-values, or equivalently, the nuclear masses of the respective nuclides, 



for reactions taking place in type I X-ray burst nucleosynthesis have been measured experimentally, 



although a number of nuclear masses need to be estimated by extrapolating experimental values [174 
For several crucial links, the reaction Q-values are currently not based on measurements (see Section 



4.4). Besides Q-values, the rates of second-step reactions in sequential two-proton capture on waiting 
point nuclei, for example, 65 As(p,7) 66 Se, 69 Br(p,7) 70 Kr and 73 Rb(p,7) 74 Sr, are expected to be important 



for nucleosynthesis considerations (Section 4.3). 

Apart from Q-values and reaction rates, decay constants for /3-decays are needed for modeling the 
nucleosynthesis in type I X-ray bursts. These have been measured with sufficient accuracy for most 
nuclides up to mass A = 100 175] (see Section 4.4). As was the case for reaction rates, the "laboratory 

For type I XRBs 



176-178 



decay constants" should generally be replaced by "stellar decay constants" 
the correction factors are presumably small for most nuclides, though detailed calculations are needed 
to better quantify these factors and examine their effects in XRB models. 



4.2 Overview and applicability of statistical model calculations 



As discussed in Section |3.3[ a typical type I XRB reaction network calculation contains thousands of 
nuclear links (reactions and weak interactions). Only a tiny fraction of the thermonuclear reaction rates 
are derived from directly measured nuclear physics input. The few direct measurements that have been 
performed, necessarily involving radioactive ion beams, focused mainly on the breakout sequences from 



the CNO mass range (Section 4.4). It is thus clear that the vast majority of reaction rates must be 
adopted from a nuclear reaction model. 

The most successful model for theoretically estimating a large number of reaction rates for nucle- 



osynthesis studies is based on the Hauser-Feshbach theory of nuclear reactions 179, 180 . It assumes 
that near the incident energy there is a large number of levels for each quantum number (spin and 
parity) in the compound nucleus through which the reaction can proceed, implying a large level density 
in the astrophysically important excitation energy range of the compound nucleus. If this condition is 
fulfilled then the Hauser-Feshbach theory requires very few ingredients for predicting cross sections or 
reaction rates. The most important ingredients are transmission coefficients and level densities. Both 
of these input quantities could be fine-tuned for a given reaction of interest, for example, by adopting 
parameters that reproduce measured cross sections. In reality, however, the number of unmeasured 
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reactions is very large and it is necessary to compute the desired cross sections and reaction rates with 
global instead of local parameters. For comparisons of Hauser-Feshbach predictions with measured cross 
sections, see Refs. fl6^|[l8T|[l82 



It is commonly assumed in nuclear astrophysics that the Hauser-Feshbach model can reliably predict 



reaction rates on average within a factor of 2 (e.g., Refs. 171,183 , but see below). Another major 
advantage of the Hauser-Feshbach model is that the effects of thermal target excitations on the reaction 
rate can be included in a straightforward manner. An extensive discussion of the Hauser-Feshbach 



model in nuclear astrophysics applications can be found in Goriely et al. |172| and Rauscher 184 . 

How well do current Hauser-Feshbach rates for charged-particle reactions reproduce the rates based 
on experiment? A preliminary answer is given in Fig. |9j showing the ratios of theoretical and experimen- 
tal rates versus both mass number (top panel) and projectile separation energy (bottom panel). The 
experimental rates, including statistically rigorous uncertainties, are adopted from a recent Monte Carlo 



evaluation 185 . Rates calculated using the code NON-SMOKER 171 are shown as filled circles, while 



those computed with the code TALYS |186| are displayed as open squares. Rates for proton-capture 
are shown in red online and those for (p,a) or (a, 7) reactions in blue online. For each reaction the 
theoretical and experimental rates are compared at the highest temperature (T match , different for each 
reaction) at which the experimental rates can be fully based on nuclear physics data. For all ratios 
shown, the "minimum temperature for applicability of the Hauser-Feshbach model" (T min , as discussed 
in Ref. |181| ), is smaller than T matc h- A few interesting observations can be made: 

• for (p/y) reactions on target nuclei in the A = 20 — 40 range, there is considerable scatter in 
the ratios of Hauser-Feshbach and experimental reaction rates (note the dashed horizontal lines, 
indicating a factor of 2 deviation), up to factors of 100 in a few cases. 

• the theoretical (p,7) rates from NON-SMOKER and TALYS also differ considerably in several 
instances, up to an order of magnitude. This is not surprising, considering that these codes 
employ different global models for the nuclear physics input (i.e., transmission coefficients and 
level densities). 

• for the five (p,a) or (a, 7) reactions shown in blue online, the agreement between Hauser-Feshbach 
and experimental results, as well as between NON-SMOKER and TALYS rates, is much better 
than for (p,7) reactions. This is especially surprising because relatively fewer levels will contribute 
to the rates of (p,a) and (0,7) reactions compared to proton capture because of angular momen- 
tum selection rules (at the same excitation energy in the compound nucleus above the projectile 
separation energy). 

• for projectile separation energies below 3 MeV, about half of the rate ratios exceed the "factor 
of 2" deviation band; i.e., there is a large risk that adopting Hauser-Feshbach results will over- 
estimate the actual rate significantly. This is not surprising because the level densities at these 
excitation energies are relatively small and the conditions for the applicability of statistical reac- 
tion models are rarely fulfilled. However, this also implies that the minimum temperature, T min , 
for applicability of the Hauser-Feshbach rates 181] ) is significantly underestimated. 



• for projectile separation energies above 5 MeV there is generally better agreement, as expected 
from level density arguments. There is some scatter, but on average the Hauser-Feshbach calcu- 
lations seem to reproduce the experimental rates within a factor of two. 

It should be clear that Fig. [9] does not represent a rigorous test of the reliability of Hauser-Feshbach 
reaction rates for type I X-ray burst nucleosynthesis simulations. First, the data set is rather small, 
especially for the (p,a) and (a, 7) reactions shown in blue online. Second, nothing can be concluded 
regarding other types of reactions, such as (a,p), or for target masses in excess of A = 40. Third, for 
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obvious reasons, the ratios presented here apply to stable target nuclei only. However, one may suspect 
that the deviations between Hauser-Feshbach and experimental rates will be larger for unstable target 
nuclei. Clearly, more work is needed to quantify the effects discussed above. 



4.3 Sensitivity studies to address nuclear physics needs 



We have mentioned the large number of nuclides and interactions needed to accurately determine the 
nucleosynthesis in type I X-ray bursts, whether within the framework of multi-zone hydrodynamic simu- 
lations or simpler one-zone models. Sufficient experimental information to calculate an experimentally- 
based thermonuclear rate over XRB temperatures only exists for a limited number of reactions. There- 
fore, most reaction rates in an XRB model depend upon theoretical calculations, as discussed in Section 
4.2 To reduce the reliance on these theoretical rates, experimentalists should measure nuclear physics 
quantities (e.g., resonance energies, resonance strengths) needed to extend existing experimentally- 
based rates to the relevant temperatures for XRBs. Furthermore, reaction rate equilibria can be better 
quantified and the input to theoretical rate calculations can be improved through measurements of, e.g., 
nuclear masses. Since nucleosynthesis in XRBs eventually extends far into the proton-rich side beyond 



the valley of stability (Section 3.3), it is clear that radioactive beams are needed to perform many of 
the desired measurements. To maximize the impact of these measurements, experiments should focus 
upon the reduction of those uncertainties in reaction rates and nuclear masses that most significantly 
influence model predictions of XRB observables. 

A useful technique to identify important nuclear interactions involves the examination of the impact 
on model predictions from the systematic variation of each interaction in a network by its uncertainty. 
Each rate may be varied independently (requiring, in theory, as many model calculations as there 
are rates in the associated network) or simultaneously through a Monte Carlo method (from which 
e.g., correlations between an enhancement factor applied to a rate X and the corresponding change in 
the yield of an isotope Y may be deduced). For theoretical rates in the network, suitable variation 
factors may be adopted to account for possible discrepancies, as discussed in Section |4.2[ These types 



of sensitivity studies have already been performed for several other astrophysical scenarios, exposing 
nuclear processes that affect nucleosynthesis in, e.g., classical nova explosions 187 188 

or massive stars 



core-collapse supernovae 190 191 



192 



AGB stars 189 



For type I XRBs, comprehensive sensitivity studies using full reaction networks coupled to multi-zone 
hydrodynamic models are not yet feasible due to computational limitations. Some progress has been 



made in this context through the use of different rate libraries |193| , variation of groups of rates 59,194 



or the variation of a few specific rates [37 92 [155 195 . For example, Cyburt et al. 193 found modest 
changes in the computed light curves (by ~ 5% in relative luminosity) and larger changes in the 
composition of burst ashes (greater than a factor of 2 in some cases) when exploring the effects of using 
two different rate libraries. Woosley et al. 
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investigated the sensitivity of their calculated light curve 
to five sets of weak interaction rates, finding large deviations in rise times, peak luminosities L pea k and 
evolution of the light curves after L pea k. Thielemann et al. 194 reported dramatic differences in light 



curves (shapes and L peak ) calculated with two networks identical except for four rates (proton-capture 
on 27 Si, 31 S, 35 Ar and 38 Ca). Fisker et al. (92] multiplied the rates of the 30 S(a,p) and 34 Ar(a,p) reactions 
by factors of 100, finding minor effects from each variation relevant to a double-peak structure in their 
light curves (see Section 3.3.2). Finally, Fisker et al. 37 155 and Davids et al. 195 , using different 
codes, varied the 15 0(a,7) rate within uncertainties to ascertain its impact on light curves. Curiously, 
contradictory results indicative of model dependencies are seen. When Fisker et al. 155 used what 
they refer to as a "lower limit" for the 15 0(a,7) rate, their model revealed non-bursting behaviour - 
more precisely, a slowly oscillating luminosity with a period of about four hours. The use of a larger 
rate produced bursts in their model, with L pea k about two orders of magnitude greater than luminosities 
they found with the lower rate. On the other hand, when Davids et al. 195 used a very similar rate 
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to Fisker et al.'s "lower limit" (within a model with similar accretion rate as well) their simulation not 
only produced bursts, but the brightest bursts; namely, a larger 15 0(a,7) rate resulted in bursts with 
~ 50% lower L pea k values. It is clear that improved tests with different models must be performed to 
clarify these discrepancies and others previously identified (e.g., the role of compositional inertia - see 
Section 3.3.2). 

Systematic sensitivity studies to guide future experimental nuclear physics measurements for XRBs 
have only been made with one-zone or post-processing models. As discussed in Section |3.3.1[ these 



are significantly less computationally expensive than multi-zone models, but, by construction, neglect 
crucial hydrodynamic effects such as convection and the finite diffusion time for energy to escape from 
the accreted envelope. Some effects of these simplified approaches on predictions of the nucleosynthesis 



and energy generation are illustrated in Figs. 10 and 11 We emphasize that the comparisons in Figs. 10 
and 11 are only strictly valid for the models discussed below. 

Fig. 10 shows the final mass-averaged yields from a multi-zone hydrodynamic model (open squares, 



burst 3 from model 1 of Jose et al. |6l|). These are compared to the final yields from a post-processing 
calculation using the thermodynamic history (see Fig. [3]) from only the hottest shell of the same hydro- 
dynamic model (filled circles). The abundances from the hydrodynamic model at T = 0.2 GK during 
the initial stages of the burst were used as initial abundances for the post-processing calculation, and 
the final yields for both calculations were extracted at T = 0.2 GK during the decline of the burst. 
Discrepancies are clear, particularly below A m 25; nonetheless, at higher masses the results for both 
calculations show similar trends, especially in the higher relative abundances at the waiting points (A 
= 60, 64, 68, 72, 76, 80). 

Fig. [TT] shows the light curve from the same hydrodynamic model mentioned above, the nuclear 



energy generation rate within the hottest shell of the same model, and the nuclear energy generation 
rate calculated during the post-processing of the thermodynamic history (see Fig. |3| from the same shell 
of the same hydrodynamic model. The energy generation maxima occur roughly where the waiting- 
point abundances change rapidly - see e.g., Ref. 151 . Although one may be able to compare selected 



characteristic features in the nuclear energy generation curves, it is evident that any careful study 
of changes to the nuclear energy generation rate (let alone a predicted light curve) owing to nuclear 
physics uncertainties can only be made with a full hydrodynamic model. The energy contributions 
from different shells as well as the time required for energy to leave the envelope clearly make non- 
negligible contributions to the light curve. As well, we stress again that convection is not included in 
these post-processing calculations. 

Nonetheless, if a post-processing sensitivity study relies upon a large number of thermodynamic 
histories, representing the available parameter spaces of temperatures and densities achieved in an 
XRB, the results serve as a useful guide for future measurements. This is especially true if variations 
in a particular reaction rate or Q-value affect predicted quantities in multiple one-zone models (where 
each uses a different thermodynamic history). Such an approach is necessary in light of the wide variety 
of observed XRB light curves, possibly indicating diverse conditions under which the thermonuclear 
runaway occurs. The usual note of caution with regard to these types of studies is worth stressing: 
since post-processing calculations only track existing thermodynamic histories, results obtained from 
variations in nuclear processes that significantly affect the energy production should be interpreted 
carefully. Indeed, a self-consistent analysis with a hydrodynamic code capable of suitably adjusting 
both the temperature and the density of the stellar envelope (in response to the change in energy 
generation) seems necessary to reliably treat such cases. 

Some limited investigations of the impact of reaction rate uncertainties on predictions of one-zone or 
post-processing models have been made previously. For example, Iliadis et al. 151 varied independently 
the rates of proton-capture reactions on 27 Si, 31 S, 35 Ar and 39 Ca and found only modest variations in 
final yields of ~ 50% or less for masses A = 40 - 48, and negligible changes to the total energy generation. 
Schatz et al. 159 and Koike et al. 140 examined the impact on yields and energy production from 
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using different reaction rate libraries. Wallace and Woosley |137| and Koike et al. 141 gauged the 
impact of the size of their networks on energy generation and the extent of the nucleosynthesis. Indeed, 
Koike et al. 141 found that XRBs may produce species up to 126 Xe, beyond the SnSbTe cycle; they 
ascribe the necessary abundance flows at the heaviest masses to their larger network compared to that 
of Ref. 153 . Clement et al. 196 reported considerable differences in light curve shapes resulting from 
calculations with different nuclear mass models (for nuclei A = 60 - 75); these are essentially effects 
due to different separation energies around waiting point nuclei, as discussed in Section |4.1| As well, 



several one-zone calculations have recently been performed to examine the impact of particular new 
experimental results (see Section 4.4). 

Post- 



144, 145 



The only full sensitivity studies for type I XRBs were reported in Parikh et al. 
processing calculations were performed using seven different thermodynamic histories, comprising mod- 
els extracted from the literature and parameterized models. These scenarios were chosen to explore 
the range of temperatures and densities achieved during XRBs and find common rates and reaction 
Q- values whose uncertainties needed to be reduced. As well, one of these histories was tracked assuming 
initial abundances of different metallicities to explore the impact of this variable. Three different ap- 
proaches were employed: first, the rates of all nuclear processes within the network were independently 
varied within uncertainties, for each scenario; second, rates for all reactions with Q-values less than 
1 MeV were recalculated assuming Q-value variations within uncertainites, and the impact of each of 
these variations was examined for each scenario; and third, rates of all nuclear processes were varied 
simultaneously using a Monte Carlo procedure, for two scenarios. An uncertainty of a factor of 10 was 
assumed for all purely theoretical rates in the network. In the independent- variation study, uncertainties 
in relatively few nuclear processes were observed to affect the final nucleosynthesis significantly (i.e., to 
change the yields of species with final mass fractions > 10~ 5 by at least a factor of two). The impact 
of variations in the 61 Ga(p,7) and 65 As(p,7) rates is illustrated in Figs. 



12 



and 13. These are two ex- 



amples of rates that did significantly modify the nucleosynthesis in most of the adopted scenarios; only 
theoretical calculations are available for both of these rates. The Q-value variation study examined the 
impact of uncertainties in nuclear masses on abundance flows about nuclear waiting points. A principal 
result of this part of the study was the need to determine experimentally the proton separation energy 
As (see Section 4.4). Figure 14 (thin grey bars) shows the dramatic effect on final yields found 
65 As) within the uncertainty estimated in Ref. 



174 . Finally, the results of the Monte 



S p of 65 

from varying S p 

Carlo and independent- variation approaches were in agreement: similar reactions were identified with 
rate uncertainties that significantly affected final XRB yields. 

Given the recent availability of a new set of hydrodynamic models [61], we have repeated the 
independent-variation sensitivity study using a thermodynamic history from one of these bursts (see 
Fig. [3]). In general, the results are in accord with those reported in Ref. |144] . Variations in some 



reaction rates 
Figs 



e.g., Ga(p,7) and 65 As(p/y)) have similar effects on extended mass regions (as in e.g. 



12 



and 



13) 



: 



In addition, variations in the rates of alpha-capture reactions on 24 Mg, 28 Si, 29 Si and 
P(p,7) 31 S reaction were found to affect the final yields of at least three isotopes around 



32 S and the 

A « 30 by at least a factor of two (similar to results seen with scenarios F08, K04-B2, K04-B4, and 



K04-B5 in Ref. [144] ). 



We note in passing that variations in the 15 0(a,7), 27 Si(p,7), 31 S(p,7) and 35 Ar(p,7) rates by exper- 
imental uncertainties did indeed affect the nuclear energy generation in the post-processing calculations 
of Ref. |144| (and in the new set of calculations described above). As such, hydrodynamic tests seem 
required to examine the impact of uncertainties in these rates. As mentioned above, several studies 
of the 15 0(a,7) rate uncertainty have already been made, but new hydrodynamic calculations of the 
impact of uncertainties in the other rates seem warranted given the results of Refs. 



as new experimental results that may affect these rates 210 



151,194 , as well 



As full sensitivity studies with multi-zone hydrodynamic models are not yet possible, and one- zone 
or post-processing calculations suffer from necessary simplifications, an obvious compromise would be 
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to first identify reactions of importance through the latter approaches and then investigate variations 



in only those selected reactions through hydrodynamic models. Such studies are in progress 1 142 , 143 



however initial results indicate that this method may not be a satisfactory substitute for a full sensitivity 
study using only hydrodynamic models. For example, Amthor et al. 142 found conflicting results 
when performing rate variations in one-zone and hydrodynamic models: variation of the 46 Cr(p,7) rate 
produced similar results with the two approaches, but differences in the results appeared when the 
42 Ti(p,7) or 49 Fe(p,7) rates were varied. As well, hydrodynamic tests with the processes identified in 
Refs. 



144 



145 1 are already underway. 



4.4 Recent experimental efforts to improve model predictions 



A number of nuclear physics experiments have been performed in recent years to better define model 
predictions of nucleosynthesis and light curves for type I XRBs. The majority of these are high- 
precision mass measurements, highlighted by the determinations of masses to experimentally constrain 
the proton separation energies of 65 As and 69 Br for the first time (—90 ± 85 keV and —785^40 keV 
respectively 197 , 199 201] , to be compared to the estimated values and uncertainties of —80 ± 300 keV 
and —450 ± 100 keV 174 ). The impact of the reduction in the uncertainty of S p ( e5 As) was explored 



through a one-zone model 1 20 1 , where the precision of light curve and nucleosynthesis predictions was 



seen to be improved. This is also seen in Fig. 14 , where the thick black error bars indicate the variation 



in final yields from post-processing calculations in which S p ( As) was varied within its experimental 
uncertainty. These results are in agreement with those presented in Ref. [20 1] : only immediately about 
A = 64 are the predictions still sensitive to the uncertainty in S p ( 65 As) (and even these yields change by 
less than a factor of two). This is a key issue when considering the required precision of a mass value to 
constrain nucleosynthesis predictions: the answer clearly depends upon which mass is being measured 
and which mass region is being examined in the resulting XRB yields. This has been addressed by 
Schatz et al. 202 in the context of changes to the effective lifetimes for 64 Ge, 68 Se and 72 Kr from 
uncertainties in nuclear masses, and in the calculations of Parikh et al. 145 where variations of a few 



specific reaction Q-values were made by then-current uncertainties and 20% of those uncertainties. For 
example, it was found that varying S l p ( 64 Ge) by ± 300 keV affected a large number of yields (thin grey 

60 keV did not change any yields by more than a 



error bars in Fig. 14) while variations by 0.2x300 
factor two. 

The expected effect of the measured S p ( 65 As) and S p ( 69 Br) values would be a reduction in the 



abundance flow beyond A fa 70 |200| relative to predictions using the estimated masses from Ref. 174 



This is supported by preliminary hydrodynamic calculations using model 1 from Ref. 61 , although the 



influence of the new measurements on the mean composition of the envelope following a burst is not 
large. Indeed, we find reductions in mass fractions of less than 30% for nuclides A > 65. Additional 
hydrodynamic models are needed to explore further the precise impact of these measurements in detail. 

Other recent relevant mass measurements have provided precise values for masses A > 80, including 
proton- rich isotopes of elements between Rb and I 154 203 208f 



Some of these studies have been 



motivated in part by exploring the still-unclear maximum endpoint of nucleosynthesis in XRBs 154 
or explaining the solar abundances of the light p-nuclei 92,94 Mo and 96 ' 98 Ru through contributions from 
XRBs (see Section 3.3.2). At lower masses, measurements have improved the S p ( 57 Cu) value 209 
and the 27 Si(p,7) rate through a new measurement of the 28 P mass [210] . From the list of masses 



desired to constrain nucleosynthesis in the post-processing calculations of Ref. 144 , 145 , experimental 
determinations of the masses of 26 P, 27 S, 43 V, 46 Mn, 47 Mn, 51 Co, 56 Cu, 62 Ge, 66 Se, 7U Kr, 84 Nb, 86 Tc, 89 Ru, 



See also databases at http://www.nuclearmasses.org/ 



http : / /research . j yu . fi/igisol/ J YFLTR AP .masses / 
http://www-win.gsi.de/aptrap/database/index.asp 
http: / /isolt rap . web . cern. ch /isoltrap/database/isodb . asp 
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90 Rh, 96 Ag, 97 Cd, "In and 103 Sn are still needed. In addition, the experimentally-known masses of 31 C1, 
45 Cr, 61 Ga, 71 Br and 83 Nb should be determined to better precision. These desired mass measurements 
are summarized in Table [U 

With regard to lifetimes of radioactive nuclei involved in XRB nucleosynthesis, the required /3 + decay 
and electron capture rates depend upon the temperature and density of the environment. While such 
"stellar" weak interaction rates have been computed for different ranges of nuclei (e.g., A = 1 
A = 45 - I 

'laboratory 



39 177 



178 ; both using the shell model) no single consistent treatment of stellar lifetimes for all 

and 



isotopes in an XRB nucleosynthesis network exists. Although the differences between 
"stellar" lifetimes are presumably small given the temperatures involved in XRBs (certainly compared 
to supernovae) the effects of these differences on model predictions need to be examined in more detail. 
We note that the vast majority of "laboratory" lifetimes for the nuclei of interest in XRBs have not 
only been measured, but indeed, are claimed with uncertainties of less than ^30% 175 . To illustrate, 
even for the species with unknown masses listed above, only the lifetime of 51 Co is unknown; lifetimes 
of the other nuclei are all known to ~ 30%. Some of these values may however be based upon only one 
measurement and so efforts to ensure the reliability of these evaluated lifetimes are welcome, especially 

±30% or 



at waiting points 



161,212 



Present experimental uncertainties in "laboratory" lifetimes (i.e., 
better) were found to affect significantly neither any XRB yields nor the nuclear energy generation in 
the post-processing studies of Ref. 



144 



214 215 



p(p, 7 ) po , 3o s(p, 7 ) m 



A few direct or indirect experiments have been performed recently to improve several reaction 
rates for XRB calculations, including the breakout reactions 14 0(a,p) |213 , 15 0(a,7) |156| . 

and 21 Na(p,7 ) |216| and reactions around the 30 S waiting point, such as 29 P(p,7) 

S(a,p) (221], 32 Cl(p, 7 ) (222), and 35 Ar(p, 7 ) [l98). 



18 Ne(a,p) 
2T7H219I, 



30 c 



5 Outlook 



Tremendous progress has been made in recent years to better understand nucleosynthesis in type I X- 
ray bursts and the role of nuclear physics uncertainties in model predictions of observables. Prospects 
for stronger observational tests of XRB nucleosynthesis predictions depend upon future X-ray missions 
and the development of improved models for the atmospheres of neutron stars. The Rossi X-ray 
Timing Explorer had been the most versatile observatory for XRBs over the past decade because of its 
ability to monitor the sky for transient outbursts from low-mass X-ray binaries and subsequently trigger 
pointed observations with its higher resolution instruments. Fortunately, other missions still in orbit 

and soon, eROSITA [225] ) 
Chandra (TOO], Swift /XRT 



have wide field cameras onboard (Fermi [223], Swift/BAT (29], MAXI (224 
while detailed pointed measurements are available with XMM-Newton |28| 
or Suzaku 
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Some recent missions such as NuSTAR 227] (launched in June 2012) and future, 
proposed projects such as the Large Observatory For X-ray Timing (LOFT |228| ) will help to advance 
our understanding of XRBs. The Large Area Detector planned for LOFT would have an effective area 
of ~12 m 2 , which is an order of magnitude larger than current spaceborne X-ray detectors (in the 2 - 
30 keV energy range). This instrument, as well as NuSTAR 229 , should be able to detect absorption 
edges in PRE bursts. The second instrument onboard LOFT, the Wide Field Monitor, would monitor 
transient sources over a large fraction of the sky. With these instruments, LOFT would be able to 
probe the spectral variability of sources down to unprecedentedly short timescales. Unfortunately, no 
currently planned mission, LOFT included, nor any mission currently in orbit has both the collecting 
area and the necessary energy resolution to measure absorption lines in neutron star atmospheres over 
reasonable timescales. 

Improved models of neutron star atmospheres are needed to aid in the interpretation of the observed 
neutron star spectrum. The nature of the ashes of thermonuclear burning during XRBs represents a 
key ingredient. Currently, the use of atmosphere models is limited to the determination of the colour 
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correction factor from assumed atmosphere abundances. This factor is then used to determine the 
actual neutron star photospheric temperature from the deduced blackbody temperature. Direct fitting 
of the observed spectra with a grid of atmosphere models could provide constraints on the abundances 
of species present in the neutron star atmosphere. This approach, combined with high signal-to-noise, 
time-resolved spectroscopy could, in principle, even probe changes in the chemical composition of the 
atmosphere during a burst. 

The coupling of state-of-the-art hydrodynamic models and large nuclear reaction networks has led to 
a better understanding of the complex interplay between nuclear interactions and X-ray burst properties 
(such as light curves, burst recurrence times or a values). It is expected that better agreement with 
observations will be achieved through the incorporation of the neutron star spin and magnetic field 
in models. Moreover, since aspects such as convective transport and burning front propagation can 
only be accurately modeled through multi-dimensional simulations, such studies must be performed 
whenever the necessary resources become available. In contrast with related astrophysical scenarios 
such as classical novae, no mass ejection results from current hydrodynamic simulations of XRBs. This 
aspect simplifies the treatment of the problem since purely Lagrangian schemes are not strictly needed to 
properly follow the expansion of the burst envelope. However, one should bear in mind that the modeling 
of both sites, novae and X-ray bursts, is extremely challenging because of the deflagrative nature of 



the burning front. Rigorous treatments will likely require the use of low Mach number codes 230 , and 
such studies are therefore encouraged. As well, the need to rely on multiple bursts to account for the 
long-term evolution of XRB sources makes multi-dimensional simulations extremely challenging, and 
it is expected that models with spherical symmetry will persist. Alternative ways to better constrain 
XRB nucleosynthesis should also be explored. This should include the detailed modeling of the process 
of photospheric radius expansion, the possibility of the ejection of nuclear-processed material through 
radiation-driven winds and the imprint of these winds on the burst spectrum. 

Several radioactive beams facilities capable of making nuclear physics measurements needed to im- 
prove estimates of reaction rates in XRBs are online or will be online shortly (e.g., TRIUMF-ISACQ 
RIKEN-RIBFQ GSI-FAIR[5 GANIL-SPIRALSQ MSU-FRIfQ. Through post-processing calcula- 



tions, nuclear physics uncertainties that may most significantly affect model predictions have been 
identified. Measurements of the masses listed in Tabled] are needed to constrain abundance flows about 
waiting-point nuclei or improve theoretical rate calculations. We stress the need for nuclear structure 
information above the proton thresholds in 62 Ge, 66 Se and 70 Kr to better determine the rates of the 
61 Ga(p,7), 65 As(p,7) and 69 Br(p,7) reactions which act upon equilibrium abundances following waiting 
points. Resolving these reaction rates should be priorities for experimentalists. In addition, measure- 
ments to further improve the rates of breakout reactions such as 14 0(a,p), 15 0(a,7) or 18 Ne(a,p); rates 
around A = 30, particularly of interest for multi-peaked structures in predicted light curves; and others, 
such as 57 Cu(p,7) or 96 Ag(p,7) would be welcome. Furthermore, since "stellar" decay constants should 
ideally be used in XRB models rather than "laboratory" decay constants, we encourage the develop- 
ment of improved, consistent treatments for calculating stellar weak interaction rates for all isotopes in 
a typical XRB network. Finally, comprehensive sensitivity studies employing hydrodynamic models are 
needed to directly gauge the impact of nuclear physics uncertainties on observable properties of type I 
X-ray bursts. 



9 http:/ /www.triumf.ca/ 
10 http:/ /www. nishina.riken.jp/RIBF/ 
11 http:/ /www.gsi.de/en/research/fair.htm 
12 http:/ /www.ganil-spiral2.eu/ 
13 http:/ /www. frib.msu.edu/ 
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Table 1: Mass measurements desired to improve calculations of nucleosynthesis in XRBs 144, 145 



Estimated masses and uncertainties from Ref. 174 are given with a # symbol; increased precision is 



required for the other, experimental masses listed. Masses required primarily to better quantify reaction 
rate equilibria at waiting point nuclei (W) or refine theoretical rate calculations (T) are indicated. 
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Figure 1: Type I X-ray burst light curves from (a) 4U 1636-536, as observed with XMM-Newton/EPIC- 



pn 28 , and (b) GS 1826-24, as observed with Swift/XRT [29]. Typical rise times and decay times vary 
between «1 - 10 s and ~10 s - several minutes, respectively. Panels (a) and (b) were generated using 
public data, available at http:/ /heasarc. nasa.gov/, Panels (1-4) show four bolometrically double-peaked 
light curves from 4U 1636-536, as observed with RXTE [To] (PCU refers to RXTE Proportional Counter 
Units). The separation of the two peaks is ~3-5 s. Credit for panels (1-4): A.L. Watts, I. Maurer, 
A&A 467, L33, 2007, reproduced with permission, ©ESO. 
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Figure 2: Evolution of the bolometric flux (top), emitting area (middle), and blackbody 
temperature T (bottom) during a type I X-ray burst of the Rapid Burster (MXB 1730-335) 
showing photospheric radius expansion. More precisely, the ordinate of the middle panel 
shows a quantity proportional to the emitting area, where Rbbody is the radius of the source 
(in km) and D is the distance to the source in units of 10 kpc. Three phases are evident: 
expansion of the photosphere, lasting « 3 seconds; contraction of the photosphere, lasting 
~ 7 seconds; and the final cooling of the envelope. Data obtained with the Swift X-ray 
Telescope [29). See Ref. 45 for more details.^ 
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Figure 3: (Top) Temperature and (middle) density versus time during the evolution of the hottest shell 
from a hydrodynamic XRB model (burst 3 of model 1 from Ref. |6l]). (Bottom) Envelope size during 
the same burst, measured relative to the core-envelope interface. Note the modest expansion of the 
envelope during the burst, achieving a maximum height of 44 m above the core-envelope interface. The 



associated light curve for this burst is shown in Fig. 11 
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Figure 4: Effect of the metallicity of the accreted material on type I X-ray bursts. (Top, filled circles) 
Mean composition of the envelope following the fourth burst computed in a model with a 1.4 M 
neutron star, accreting solar-like material (Z=0.02) at a rate of 1.75 x 1CT 9 M yr -1 [6l]; (bottom, 
solid line) the light curve associated with this burst. (Top, open squares) Mean composition of the 
envelope following the fifth burst computed in a similar model as above, but accreting material of lower 
metallicity (Z=10 -3 ) [61] ; (bottom, dashed line) the light curve associated with this burst. Accreted 
material of lower metal content drives bursts characterized by lower peak luminosities and longer decline 
times. In turn, the longer duration of the burst allows a significant extension of the nuclear activity 
towards the SnSbTe-mass region, which is clearly reflected in the mean composition of the envelope at 
the end of the burst. 



36 





72 Sr 


"Sr 




,2 Rb 






70 Ki 


' 1 




"Br 


"Br 






"Rb >»Rb l0 Rb 



■'Mi »Mi >»Mi »Mi > 2 M, 



13 Ag 



"Rl « 2 Rt > 3 Rb i4 Rl 



"Ag 



«Nb ,7 Nb "Nb >»Nb l0 Nb "Nl > 2 N1 



»«Ag 



"Ag 



"Ag 



"Ag 



Sd I08 Sd I09 Sd 



101 Ag 102 Ag 103 Ag 104 Ag 05 Ag !06 Ag 



PC l03 Pc 04 Pd 



M< "Ml MMij 



["Zr|"Zr| 



Model 114NSz 




Log (Mass Fractions) 
-1 <log X<0 
■2 <log X<-1 
-3 < log X < -2 
-4 < log X < -3 



| ]1 Ti| ] »Ti 



pMg "Mj "Mg 



|"Ne|' 7 



20 Na 2 'Na !2 Na 



, 5 . 



"Mi 



'«J!i "Ca 




Ti 4I Ti 4, Ti 50 Ti 



T peak =1.16xl0 9 K 
p= 1.59x10 s gem 3 



E nnB =8.8xl0 16 eiggr 1 s 
X= 0.037 
Y=0.21 
X[ 64 Ge]=0.12 



Figure 5: Main nuclear activity during a XRB at T pea k, shown in terms of mass fractions for the most 
abundant species (X > 10~ 4 ) within the ignition shell (~ 5.6 m above the core-envelope interface). 
These results were found for the fourth burst computed in a model with a 1.4 M neutron star, accreting 
solar-like material (Z=0.02) at a rate of 1.75 x 10 -9 M yr" 
synthesized and the modest extension of the activity up to A : 
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Figure 6: Same as Fig. |5| but for accretion of low-metallicity (Z=10 ) material. As in Fig. |5j the most 
abundant species at T pea k is 64 Ge; however, the nuclear activity here extends up to A = 107. 
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Figure 7: Same as Fig. [5j but the nuclear activity is shown in terms of the dominant reaction fluxes 
[F > 1CT 3 ). For clarity, the numerous (p, 7)-(7,p) pairs at equilibrium are not displayed. Note 
that at Tp ea fc, the nuclear flow is dominated by proton-capture reactions (exceptions being 16 0(a,7), 
17 F(a,p), 18 Ne(a,p), and 22 Mg(a, p)). A few photodisintegration reactions (not in equilibrium with the 
corresponding (p,7) processes), and a number of B and /3-delayed decays are also shown. The main 
nuclear activity proceeds far from the valley of stability above A = 45. 
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Figure 8: Same as Fig. |6j but the nuclear activity is shown in terms of the dominant reaction fluxes. 
For clarity, the numerous (p/y)-(7,p) pairs at equilibrium are not displayed. The reduced metallicity 
of the accreted matter results in a net reduction of the reaction fluxes, with the most abundant fluxes 
being three orders of magnitude less frequent than those reported for the solar metallicity model in Fig. 
[7J Note the significant number of /3-decays presented here (compare Fig. [7]) and the extension of the 
activity up to the SnSbTe-mass region. ^ 
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Figure 9: Comparison of theoretical (Hauser-Feshbach) and experimental "laboratory" reaction rates. 

The theoretical rates are obtained using the codes NON-SMOKER (filled circles |171| ) and TALYS 

(open squares; courtesy of S. Goriely (2012)), while the experimental results (including uncertainties) 

are adopted from a recent evaluation |185| . Ratios for (p/y) reactions are shown in red online, while 

those for (p,a) or (0,7) reactions are displayed in blue online. The rates are compared at the highest 

temperature (different for each reaction) at which the experimental reaction rates can be fully based 

on nuclear physics data (i.e., at T matc h, in the parlance of Ref. 185| ). (Top) Reaction rate ratios versus 

mass number of the target nucleus. (Bottom) Reaction rate ratio versus projectile separation energy 

(equal to the Q- value for particle capture). The horizontal dashed lines represent ratios of a factor of 2 

up or down. Note the large number of values, for both computer codes, that fall outside the commonly 
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Figure 10: Comparison of final yields from a hydrodynamic XRB model (open squares, burst 3 from 
model 1 of Ref. [6l]) and from a post-processing calculation (filled circles) using the hottest shell of the 
same model. 
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Figure 11: (Top) Light curve from a hydrodynamic XRB model (burst 3 from model 1 of Ref. |61|). 
(Middle) Nuclear energy generation rate from the hottest shell of the same hydrodynamic model. (Bot- 
tom) Nuclear energy generation rate from a post-processing calculation using temperature and density 
profiles (Fig. |3| from the same shell of the same hydrodynamic model. These plots show that the rise in 
luminosity is due to the nuclear energy released within the ignition shell, deep inside the envelope, which 
is later convectively transported towards the surface. The huge spike in the nuclear energy generation 
from the hydrodynamic model (middle panel) is due to the sudden injection of protons in the vicinity 
of the ignition shell following the onset of convection. Such an effect cannot be handled within one- zone 
post-processing calculations (bottom panel). 
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Figure 12: The effect on final XRB yields from varying the thermonuclear rate of the 61 Ga(p,7) 62 Ge 
reaction by a factor of 10 (up and down). Only theoretical calculations exist for this rate. The horizontal 
bars indicate mass fractions using the standard rate, and the vertical error bars indicate changes to 
these mass fractions due to the rate variations. From post-processing calculations with model K04-B6 
in Ref. pi] . 
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Figure 13: Same as Fig. 12 but from varying the thermonuclear rate of the As(p,7) Se reaction by a 



factor of 10 (up and down). From post-processing calculations with model K04 in Ref. 144 
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Figure 14: The effect on final XRB yields from varying the reaction Q- value of the 64 Ge(p,7) 65 As reaction 
within uncertainties. The horizontal bars indicate mass fractions using either Qameo3 = —80 ±300 keV 
or Qexp — —90 ±85 keV (the former uses the estimated value and uncertainty from Ref. 174 , while the 



latter uses experimental values for all masses 197,201 . The vertical error bars indicate the changes 



to these mass fractions due to the uncertainty in Qameo3 (thin grey bars) or the uncertainty in Q 



exp 



(thick black bars). From post-processing calculations with model K04 |144 , 145 
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